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Abstract 

In this paper, we extend the study of bivariate generahsed beta type I and II 
distributions to the matrix variate case. 

1 Introduction 

Matrix variate beta type I and II distributions hav e been studied b y different author s utilis- 
i ng di v erse approa ches, see'Olkin and Rubin' ('l 964l). Khatri ( 1970VMuirhcad' ('l982)_, Cadet 



(| 19961 ) , iGupta and Nagar. (2000 ) , Di'az-Garci'a and Gutierrez-Jaimez. (2007. .2006.. JOOl), amon. 



many others. These distributions play a very important role in various approaches to proving 
hypotheses in the context of multivariate analysis, including canonical correlation analysis, 
the general linear hy poth esis in MANOVA and multipl e matrix variate correlation analysis, 
see iMuirheadI ( 19821 ) and Srivastava and Khatril ( 1979f ) . All these techniques are based on 



the hypothesis that some matrices A and B are independent with Wishart distributions. In 
the present paper, these results are generalised, assuming A and B to have a matrix variate 
gamma distributions. 

A : m X TO is said to have a matrix variate gamma distribution with parameters a and 
m X m positive definite matrix 0, this fact being denoted as A ~ Qmio-, ®), if its density 
function is 

477^|A|"-(™+i)/2g^^(_0-iA)(dA), A>0, (1) 



r™[a]|e|' 



where 



(dA) = /\day, (2) 

i<j 
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see Muirhead ( 1982 . pp. 57 and 61) and Gupta and Nagai I 200(]( ): where T„i[a\ denotes the 
multivariate gamma function and is defined as 



r,„H = / ctr(-v)|vr-('"+i)/2(dV), 

"'V>0 



Re(a) > (m - l)/2 and etr(-) = exp(tr(-)). 

As well as the classi ficati on of the beta distribution, as beta type I and type II (see 
lOupta and NagaJ ( 20001 ) and lSrivastava and Khatri (ll979D 'l. two definitions have been pro- 
posed for each of these distributions. Let us focus initially on the beta type I distribution; if 
A and B have a matrix variate gamma distribution, i.e. A ~ 5m(a, I,„) and B ^ Gmib, Im) 
independently, then the beta matrix U can be defined as 



U = 



(A + B)-i/2A((A + B)-i/2)', 
Ai/2(A + B)-i(Ai/2)', 



Definition 1 or, 
Definition 2, 



(3) 



where C^/'^(C^/'^)' = C is a reasonab le n qnsingular f actori zation of C, see lGupta and Nagar 
( 200(]| ). Srivastava and Khatri ( 19791 ) and Muirhead ( 1982[ ). It is readily apparent that under 
definitions 1 and 2 its density function is denoted as i3Jm(U; a, b) and given by 



Pm[a,b] 



|U|(a-(™+l)/2|I^ _U| 



h-(m+l)/2 



(dU), < U < I„, 



(4) 



this being denoted as U ~ .B/„i(a, 6), with Re(a) > (m — l)/2 and Re(6) > (m — l)/2; where 
(3„i[a,b] denotes the multivariate beta function defined by 

I3m[b,a] = / |Sr-(™+i)/2|i„^_S|''-("+i)/2(dS) 

"'0<S<I„ 

= f |R|"-(™+l)/2|I^^+R|-(fc+fc)(rfR) 
JR>0 

_ rm[a]rm[^] 

Tm[a + b] 

A similar situation arises with the beta type II distribution, with which we have the 
following three definitions: 



B-i/2A(B-i/2)', Definition 1, 
Ai/2B-i(Ai/2)', Definition 2, 



(5) 



with the distribution being denoted as F ~ BIIm{a, b). In this case under definition 1 and 
2, the density function of F is denoted BIIm{F; a,b) by and defined as 



1 



/3,„[a, 6] 



|p|a-(™+l)/2|j^^p|-(a+6)(^p)^ F > 0. 



(6) 



(dF) is given in analogous form to 

Some of these generalisations from a univariate beta distribution to the matrix variate 
case are inappropriate because, in some a pplications, the rese a rcher is interested in a vector 
variate, not in a symmetric matrix, see iLibbv and Novick I (|l982l ). In other words, the 
researcher is interested in a vector, say, X = . . . , a;™)', su ch that Xj has a rnargina l 
beta type I or II distribution for alH = 1, . . . , m. In this respect, Libbv and Novick ( 1982h 
and I Chen and Novick I (jl984l ) proposed a multivariate version of the beta type I an d II 
distributions. Let u s consider the following bi variate version, see lOlkin and Liul ( 20031 ) and 
Nagar et all (|2008l ). 
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Let Xq, Xi, X2 be distributed as independent gamma random variates with parameters 
a = Oo, fli, 02, respectively and = 1 in the three cases, and define 

Xi X2 



Xi + Xq ' 



U2 = 



X2 + Xq 



Clearly, Ui and U2 each have a beta type I distribution, Ui ^ BIi{ai, uq) and U2 ~ 
BIi{a2,ao), over < ui,U2 < 1. However, they are correlated such that {111,1/2)' has a 
bivariate generalised beta type I distribution over < mi, M2 < 1. The joint density function 
of Ui and U2 is 



^(1 - - ^2)"^+""-^ 



/?! [ai, 02, ao](l ~ U1U2] 



ai+a2+ao 



< Ui, U2 < 1 



where 



rm[a]r„i[6]rm [c] 



r,n[a + b + c] 

A similar result is obtained in the case of beta type II. Here define 

^ X-\ „ Xo 

Fi = — - F2 ^ — 

Once again it is evident that Fi and F2 each have a beta type II distribution, Fi ~ 
i3//i(ai,ao) and F2 ^ BIIi{a2,ao), over /i,/2 > 0. As in the beta type I case, they 
are correlated such that {Fi,F2)' has a bivariate generalised beta type II distribution over 
/17 /2 > 0. The joint density function of Fi and F2 is 



02-1 

2 



^i*[ai,a2,ao](l + /i + /2)'^ 



1+02+00 



/l,/2>0 



Some appli cations to utility modelling and Bayesian analysis are presented in lLibbv and Novick 



(|l982r ) and IChen and Novick I (|l984l ) , respectively. Properties such as the moments 
conditional distribution, t he distributions of the pro d uct u-\U2, and the mi Iu2 a nd u-\ |{u^ 



U2) qu otien ts are studied in Libbv and Novick! ( 1982I ). Chen and Novick (1984), Olkin and Liu 



(|2003l ) and lNagar et al. I (|2008l ) 



In the present paper, we extend the bivariate generalised beta type I and II distributions 
to the matrix variate case, see Section |3] and [H These distributions are termed as bimatrix 
variate generalised beta type I and II distributions. In Section [5l some properties of these 
distributions are studied. 



2 Preliminary results 

In this section, some results for the hypergeometric function with a matrix argument are 
shown. 

Definition 2.1. The hypergeometric functions of a matrix argument are given by 

where denotes the summation over all the partitions r — ti > ■■■ > 

tm > 0, of t, Cr(X) is the zonal polynomial of X corresponding to r and the generalised 
hypergeometric coefficient (a)T is given by 

m 

(a), = [](a-(i-l)/2),,, 

i=l 
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where (a)t — a{a + l)(a + 2) • • • (a + 1 — 1), (a)o — 1- Here X, the argument of the function, 
is a complex symmetric m x m and the parameters a^, &j are arbitrary complex numbers. 

Some other cha racteristics of the parameters and bj and the convergence of ([7]) appear 
in lMuirheadI (|l982l . p. 258). 



A special case of (O is 

ii^o(«;X) = (11X11 <1) 

t=0 T 

where ||X|| denotes the maximum of the absolute values of the eigenvalues of X. 

An interesting relation is derived from the following Lemma [2.H which gives an induction 
method for constructing hypergeometric functions, i.e. integration involving pFq leads to 
the new hypergeometric function p+iF^+i. A motivation for t he general recursion comes 
from the following well-known expressions, see iMuirheadI (|l982l Theorem 7.4.2, p. 264): 



iFi(a;c;X) = 



1 



f3m[a, c~ a] 

oi^o(XY)|Y|'^-(™+i)/2|i _ Y|'=~'^-(™+i)/2(rfY). 



0<Y<I„ 



and 



2Fi(a, ai;c; X) 



P,n[a, c - a] 

/ ii^o(ai;XY)|Y|"-("+i)/2|I- Y|=-'^-('"+i)/2(dY). 

JO<Y<I™, 



(8) 



Thus we have 



Lemma 2.1. Let X < I, Rc(a) > (m-l)/2, Re(c) > (rn-l)/2 andRe{c~a) > (m-l)/2. 
Then 



p+i-Pg+i (flj Q-ii ■ ■ • J Op! c, &i, . . . , &q; X) 



1 



(3.rn[a,c-a] 

pFgiai • • ■ Op; 6i • • • XY)|Yp-("+i)/2|I - Yf-'^-^'^+^^/^dY). 

/0<Y<I„ 

Proof. First, we apply an expansion in terms of zonal polynomials 

F (n n -h h ■ XY-) - V V • • • gr(XY) 

pF,(ai, . . . , ap, 6i, . . . , 6„ XY) - ^ 2^ 

Then, after integrating term by term, see iMuir head (Il98l Theorem 7.2.10, p. 254), we have 
that 

/ pFqiai ■ • ■ flp; 6i ■ • ■ 6,; XY)|Y|"-('"+i)/2|I - Y|"~"-(™+i'/2(rfY) 

JO<Y<I„ 

= V V / |Yr-(™+i)/2|i_Yr-''-('"+i)/2a(XY)(dY) 

t^O T " ' ("9)-^ "'0<Y<I„ 

= Prn[a, c- a] p+iFq+i{a, ai • ■ • Op; c, 5i • • • X), 
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and the required result follows. □ 
The use of zonal polynomials and the hypergeometric function with a matrix argu- 
ment has only recen t ly bee n exten ded; to a lar ge extent this is a derived from the work of 
Koev and Edelman ( 20061 ). who in KoevI ( 2004 ). provided a program in MatLab with a very 



efficient algorithm for calculating of Jack polynomials (in particular zonal polynomials) and 
the hypergeometric function with a matrix argument. 

3 Bimatrix variate generalised beta type I distribution 

Let A, B and C be independent, where A ^ tJm(a,Im), B ^ Qm{b,lm) and C ^ Gm{c,lm) 
with Re(a) > (to - l)/2, Re(6) > (to - l)/2 and Re(c) > (m - l)/2 and let us define 

Ui = (A + C)-i/2a(a + C)-i/2 U2 = (B + C)-i/2b(b^C)-i/2 

Of course, Ui ~ BIm{a,,c) and U2 ~ BIm{b,c). However, they are correlated such that the 

distribution of U = (Ui:U2)' G sjjsmxm termed a bimatrix variate generalised beta 

type I distribution, denoted as U ~ BQBl2my.m{p,, b, c). 

Theorem 3.1. Assume that U ^ BGBl2mxm{o-,b, c) . Then its density function is 

|U^|a-(r»+l)/2|u^|b-(r»+l)/2|l^^_U^|b+c-(,n+l)/2|l^^^_U^|a+c-(„,+ l)/2 ^^^^^ 

/3*„[a,6,c]|I„-UiU2h+''+^ ^^""^ 
< Ui <lm, < U2 < Im, where the measure 

(dU) = (dUi) A (dU2). 
andRe{a) > (to, -l)/2, Re(6) > (to-1)/2 andRe{c) > (to-1)/2. 
Proof. The joint density of A, B and C is 

I A |a-(m+l)/2|g|6-(m+l)/2|(-;|c-(m+l)/2 

^ W ' etr(-(A + B + C)){dA){dB){dC) 

By effecting the change of variable ([9]) , then 

(dA)(dB)(dC) = |C|™+i|I„, -Uir("+i'|I„, -U2|-('"+i)(dUi)(dU2)(dC). 
The joint density of Ui, U2 and C is 



|U^|a-(™+l)/2|U2|''-(™ + l)/2 



|Q|a+fc+c-(m+l)/2 



r„[a]r„[6]r,„[c]|i,„-Uih+("+i)/2|i„-U2r+("+i)/2' 

X etr [-(I™ - U2)-i(I™ - UiU2)(I™ - Ui)-iC] (dC)(dUi)(dU2). 
Integrating with respect to C using 

f |c|a+6+c-(™+l)/2 _ U2)-1(I™ - UiU2)(I,„ - Ui)-lC] (dC) 

Jc>o 

= r[a + 6 + c]J 



/c>o 

ii™-uir+''+^ii„-u2r+''+'= 



|i„,-UiU2h+''+- 

(from ((D)) gives the stated marginal density function for (Ui:U2)'. □ 



As in the bivariate case ( Olkin and Liu . 20031 ). the joint density PTI| can be represented 
as a mixture. Let us first note that 



llrn - U1U2 



-{a+b+c) _ 



iFo(a + 6 + c;UiU2) 



t = T 



a(UiU2) 

tl 



By substituting in (jTU]) we obtain that the joint density function of (Ui:U2)' is 

M-^M-^ {a + b+ c)r , |a-(m+l)/2|Tj |b--(m+l)/2|T _-rT |6+c-(m+l)/2 r' (Tl 11 \ 

0*^[aAc] i^^^i i^^^i qL^^Jj2r+-('"+^)/^ ^-^^,^^^\ 



Moreover 

(a + 6 + c)rr„i[6 + c]r„i[a + c 



EE' 



4 Bimatrix variate generalised beta type II distribution 

Let A, B and C be independent, where A fj,„(a,l„i), B ~ 5m(&,Im) and C ~ Gmicjlm) 
with Re(a) > (to - l)/2, Re(6) > (to - l)/2 and Re(c) > (m - l)/2 and let us define 

Fi = C-^^^AC-^^^ and F2 = C-^/^BC-^/^ (11) 
Clearly, Fi ~ BIIm{a,c) and F2 ~ BIIm{b, c). But they are correlated and so the distri- 
bution of F = (Fi:F2)' G S)ff2mxm ^^^^ termed a bimatrix variate generalised beta type II 
distribution, which is denoted as F BGBIl2mxm{o-, b, c). 

Theorem 4.1. Assume that F ^ BGBl2mxm{o-, b, c). Then its density function is 

IP |a-(m+l)/2|p |b-(m+l)/2 

1^1 LliJ fdFl ('121 

Fi > 0, F2 > 0, where the measure 

(dF) = (dFi) A (dFa). 
a7idRe(a) > (m - l)/2, Re(6) > (to - l)/2 and Re(c) > (to - l)/2. 

Proof. As an alternative to proceeding as in Theor em 13.11 Let us recall t hat if U ~ 
BIm.(a,b), then (Im. — \J)~^ — Im ~ B llm(a^b), see ISrivastava and Khatril (|l979h and 
Diaz-Garcia and Gutierrez- Jaimez (|2007f ). Then 



F = 



Fi 
F2 



(Im ~ Ul) ^ —Iji 
(Iiri — U2) — In 



with the Jacobian given by 

(dUi)(dU2) 

Also, note that 

Im — (Im + Fi) ^ 
Im — (Im + F2) 



|I„+Fi|-("+i)|I„ 



-(m+l) 



(dFi)((iF2). 



= (I™ + Fi)-i((Im + Fi) - I„0 = (I„ + Fi)-iFi 
= (I,„ + F2)-iF2, 



6 



(dW). 



Then the joint density of (Fi:F2)' is 

P*r^[a,b,c]\lm - (I™ + Fi)-iFiF2(I™ + F2)-i|"+^+^ 
The desired resuhs is foUow noting that 

|I„ + Fi|~^|I.,„ + F2|~^ ^ 1 

|I„-(I™ + Fi)-iFiF2(I„ + F2)-i| . 

Other properties of the distribution BQBIl2mxm{(i, b, c) can be found in a similar way. 

5 Properties 

In this section we calculate the moments i?(|Ui|''|U2|'') and the distributions of the product 
Z = U2/^UiU2^^ and the inverse (Uj^^iUa"^). 

Theorem 5.1. Assume that (Ui:U2) ^ BQBl2mxm{a,b,c) then 

_ f3m [a + r,b + c](3m [b + s,a + c] 



E{\Ui\'-\^2 



X ■iF2{a -\- r,b + s,a + b + c;a + b + c + r,a + b ^- c + s] !„), 



with Re(6 + r) > (m - l)/2, and Re(a + c) > (m - l)/2. 
Proo/. 



X 



Pm[tt!0, Cj Jo<Ui<I„ 

[ |U2|''+'^-('"+i)/2|I,„ - U2|°+'=-(™+i)/2|i^^ „ UiU2|-('^+''+=H'^U2)(dUi). 

"'0<U2<I^ 



/0<U2<I 

From (IHl) we have 



Pm [g + C, , b + s] /■ . |a+r-(m+l)/2|T _tt |6+c-(m+l)/2 



ig(iuiriu2n = ^'"^"V';Y^^ / |Uir+-("+i)/2|i„-ui 

X 2Fi{b + s,a + b + c;a + b + c + s;Vi){d\5i). 
Then, integrating using Lemma |2. li the desired result is obtained. □ 



Theorem 5.2. Assume that (Ui:U2) ^ BQBl2mxmio-,b, c). Then the density function of 
Z = vl^^lJiUl^^ is 

f3,n[a + c,b + c] |Z|'^-("+i)/2|I„ - z|--("+i)/2 



/3;„[a,6,c] 

and 



-2Fi(a + c, a + c; a + 6 + 2c; I,„ - Z)(rfZ) 



£;( Z ) = — — — 3^2(0, a + c, a + c;a + c + r, a + 6 + 2c;I„) 

with < Re(Z) < I„, Re(a + 6) > (m - l)/2 and Re{b + c) > (m - l)/2. 
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Proof. Consider the transformation Z = Ug^^UiUj^^, U2 = U2, with 

(dUi)(dU2) = \\J2\-^"'+'^^\dZ){dV2). 

Then the joint density of Z and U2 is 

|2|a-(m+l)/2|u^|-(a+c)|j^ _ jj^^a+c-{7n+l) /2 ^^j^ _ 2|b+c-(m+l)/2 



/3;„[a,6,c] |I™-Z|« 



+b+c 



(dZ)(dU2), 



< U2 < Z < I„. 

Now let us make the change of variable W = (I„j — Zi)^^/^{Ijn — U2 )(!.,„ — Z)^^^^ and 
Z = Z. Noting that 

(dZ)(dU2) = - Z|("+i)/2(dZ)(dW). 
Then the joint density of Z and W is 

" /3*^[a,b,c] - - Z)Wh" (dW)(dZ), 

0<Z<I™, 0<W< 

Integrating with respect to W using ([8]) with Re(Z) < Re(a + c) > (m — l)/2, and 
Re(6 + c) > (m — l)/2, we then obtain the stated marginal density function for Z. 

Now in order to find the expression for i?(|Z|''), let us first observe that by p = 2 and 
g = 1 in Lemma |2. II we obtain 



3^2(0, ai, 02; c, 61; X) 



X 



(3jn[a,c-a] 

[ 2^1(0102; 5i; XY)|Yr-(™+i)/2|i^ _Yr-''-(™+i)/2(dY). 

JO<Y<I,„ 



/0<Y<I 

Now, by making the transformation W = I„i — Y, with (dY) = (rfW). Then 

3-^2(0, ai,a2;c, 61; X) = 7 

Pm[a,c- a\ 

X / 2i^i(aia2;6i;X(I™-W))|I™-Wr-(-+i)/2|w|'=-~("+i)/2(dW). (13) 

Jo<w<i™ 

Then the expression for i?(|Z|'') follows immediately from the density function of Z and 
(HSl). □ 

Theorem 5.3. Assume that (Ui:U2)' ^ BQBl27n^m{0'ib,c). Then the density Junction of 
V=(Vi;V2)' = (Ur'iU2-i)' IS 

/3;;,[a,6,c]|I,„-(ViV2)-ih+^+^ ^ ^'^^ 

< Vi < Im, < V2 < Im; where the measure 

(dV) = (dVi) A (rfV2). 

anrf Re(a) > (m - I)/2, Re(6) > (rn - I)/2 a^d Re(c) > (m - I)/2. 

Proof. The desired result follows immediately by making the change of variable (Vi : V2)' = 

(U]^^:U2"^)' in the joint density of (Ui:U2)', taking into account that {d\5i){d'\J2) = 
|Vi|-("-i)|V2|-('"-i)(dVi)(dV2). □ 
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6 Conclusions 



With the algorithms proposed by iKoev and Edelman I (|2006l) for the calculation of Jack 
polynomials and hyper geometric fun ctions with matrix arguments, together with the Mat- 
Lab implementation by KoevI ( 2004 ). it is now is feasible to evaluate in a very efficient way 
expressions such as density functions and moments, as shown in the proceding sections, as 
well as highly complex expressions. The theory developed in this paper has not yet been 
applied, with the exception of the bivariate case. However, its potential role is apparent, 
for example, in multidimensional scaling (MDS) in the following context. Bimatrix variate 
generalised beta distributions can be used as distributions of the matrices of similarities (or 
dissimilarities) for an individual when these matrices of similarities have been obtained at 
two different times. Given the genesis of bimatrix variate generalised beta distributions, 
such distributions may allow us, in some sense, to model the learning problem. Statisti- 
cal approaches to MDS have been studied assumin g indep endence between t imes (without 
learning) in the univariate case bv iRamsav ( 1982| ) and bv lVera et al. ] llooi). 

Another potential use appears in the context of shape theory, specifically in the approach 
kn own as affine shape or con figuration densities. This approach is currently being studied, 



see 



Caro-Lopera et ail ( 2008 ). It was first proposed bv lGoodall and Mardial ( 19931 ). and at 



present only the Cauchy configuration dens ities, have been explored. T he use of beta con- 
figuration densities, only was proposed by iGoodall and Mardial (|l993l ). but no additional 
information was provided. Perhaps bimatrix variate generalised beta distributions can be 
obtained as configuration densities if we assume that two figures (images) of a single indi- 
vidual or a single object, obtained at two different times, are not independent. Currently, 
both applications are under study. 
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